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ABSTRACT 

The  physical  assumptions  and  mathematical  framework  are 
given  of  a  simplified,  2-dimensional  model  of  a  hydrodynamic 
disturbance  due  to  passage  of  a  submerged  submarine  in  an 
ocean  characterized  by  a  thin  thermocline.  Several  exploratory 
applications  of  the  model  are  described,  including  estimation 
of:  currents  at  the  surface  due  to  i)  the  complex  disturbances 

made  by  a  submarine  running  slightly  "heavy"  above  a  thin 
thermocline,  and  ii)  to  an  oscillating  and  eventually  collapsing 
wake  just  below  a  thin  thermocline,  and  of  the  energetic 
efficiency  of  production  of  the  lowest  mode  internal  wave  of  a 
thin  thermocline  due  to  wake  collapse. 
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NOTATION 

amplitude  of  hull  disturbance  -  A  sin  yt 
hull  radius^  or  initial  wake  radius 
radius  of  equivalent  fully  mixed  wake 
thermocline  thickness 

x/t  =  const,  locus  of  stationary  phase 
bottom  depth 
delta  function 

—  ~  P1  relative  density  jump  across  thermocline 
ps 

source  depth  below  thermocline 
acceleration  of  gravity 

frequency  factor  in  representation  of  hull  disturbance 
as  A  sin  yt 

Heaviside  unit  function 
thermocline  depth 
£s  -  p(s) 
hull  length 

strength  of  moving  source  representing  final  collapse 
stage  =  M/t 

local  Vaisala  frequency 

frequency  factor  in  time-dependent  quadrupole  strength 
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cp('r)  time-dependent  source  function 
px  density  of  upper  layer 

p2  density  of  lower  layer 

s  wave  number,  variable  of  integration 

SQ  root  of  da/ds  =  jB 

a  factor  in  integral  defining  velocity  due  to  unit  move. 

Equations  [l]  and  [2] 

a  value  of  o(s)  for  s  =  s 

o  v  '  o 

t  time 

T  interfacial  Vaisala  period  =  2tt/N 

Ti  =  2ir/y  duration  of  source  representing  hull  disturbance. 
Equation  [4] 

t  source  time  parameter 

9  factor  in  integral  defining  velocity  response  to  unit 

source.  Equations  [l]  and  [2] 

u  (x,t$  ?,t)  velocity  at  surface  due  to  unit  source  at 

(5*  ~f^T) 

U  (x,  h,  t)  velocity  at  surface,  at  time  t  due  to  time 

dependent  source  function  to(t) 

Ux  integrated  velocity  at  surface  due  to  hull.  Equation  [41 

U2  integrated  velocity  at  surface  due  to  wake  collapse. 

Equations  [7I  and  [8] 

U2  secular  (non-propagating)  part  of  U2 ,  Equation  f 7 1 
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U2  wavelike  (propagating)  part  of  U2 ,  Equation  [8] 

d3  integrated. velocity  at  surface  due  to  final  phase  of 

\ 

wake  collapse 

U3  secular  part  of  U3 ,  Equation  T9] 

U3  wavelike  part  of  U3 ,  Equation  [10] 

V  hull  velocity 

¥  source  speed  for  final  stage  of  collapse 

?  x  coordinate  of  source 

X  coordinate  along  axis  of  track  of  submarine 
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INTRODUCTION 

For  several  years,,  use  has  been  made  at  HYDRONAUTICS, 
Incorporated  of  a  simplified  model  of  hydrodynamic  disturbances 
simulating  those  that  would  be  produced  by  passage  of  a  sub¬ 
merged  submarine  in  ocean  conditions  characterized  by  a  thin 
thermocline  (more  accurately,  pycnocline).  A  brief  description 
of  some  early  results  obtained  with  this  model  were  given  by 
M.  P.  Tulin  in  Reference  1.  Recently,  more  sophisticated  models 
have  been  introduced  with  which  analytical  estimates  and  numer¬ 
ical  calculations  habe  been  made  of  hydrodynamic  effects  due  to 
internal  waves  excited  by  translational  movement  of  the  submerged 
hull  and  collapse  of  the  submarine  wake.  Recognizing  that  these 
models  can  give  more  accurate  results,  we  have  continued  to  make 
use  of  our  more  simplified  model  for  exploratory  purposes,  to 
obtain  orders  of  magnitude  of  effects  such  as  surface  currents 
and  strains  due  to  sources  to  which  less  attention  has  been  paid, 
such  as;  wakes  with  net  vertical  momentum;  and  to  combinations 
of  this  with  the  complex  array  of  sources  simulating  the  sequence 
of  hydrodynamic  disturbance  associated  with  passage  of  a  sub¬ 
merged  submarine.  It  is  hoped  that  the  results  will  be  suffic¬ 
iently  suggestive  to  stimulate  further  work  along  these  lines 
with  more  accurate  models. 

Given  below  are:  1)  a  description  of  the  physical 
assumptions  of  the  model;  2)  a  brief  outline  of  its  mathematical 
development;  3)  comparison  of  some  of  its  results  with  experi¬ 
mental  data  and  with  other  theoretical  results  for  the  case  of 
a  thin  thermocline;  4)  results  of  applications  of  the  model  to 
the  case  of  a  submarine  moving  "heavy"  in  the  mixed  region  above 
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the  thermocline,  and  moving  "light"  under  the  thermocline;  and 
to  the  (energetic)  efficiency  of  wake  collapse  as  a  source  of 
lowest  mode  internal  waves;  and  5)  a  conclusion,  including  some 
remarks  regarding  limitations  of  the  model,  and  in  regard  to 
dependence  of  its  results  on  operational  and  environmental 
parameters . 

1.  PHYSICAL  ASSUMPTIONS 

The  model  embodies  several  simplifications.  First,  it  is 
two-dimensional,  and  may  be  valid  to  the  extent  that  a  "slender 
body"  approximation  holds,  i.e.,  the  flow  parameters  change  only 
gradually  along  the  direction  of  hull  motion  and  wake  length, 
and  relatively  rapidly  in  the  plane  perpendicular  to  this  direct¬ 
ion.  For  wake  phenomena  at  all  but  low  speeds,  such  an  approxi¬ 
mation  may  be  valid.  For  disturbances  due  to  the  hull  itself, 
the  approximation  may  be  less  appropriate. 

Point  source  hydrodynamic  disturbances  are  assumed,  with 
strengths  varying  with  time  corresponding,  approximately,  to 
phenomenological  features  of  the  actual  sources.  Thus  the  dis¬ 
turbance  due  to  a  [  cylindrically  symmetric)  hull  of  radius  a 

o 

and  length  L  moving  with  velocity  Vq  is  represented  by  a  two- 

dimensional  monopole  proportional  to  sin  yt  where  y  =  2wV  /L; 

o 

the  early  stages  of  collapse  of  a  wake  of  radius  a  in  a  density 

gradient  with  Vaisala  frequency  N  is  represented  by  a  point 

quadrupole  of  strength  proportional  to  te-vt  where  \>  „  K  is  the 
frequency  corresponding  to  the  maximum  of  the  collapse  spectrum 
(see  Wu,  Reference  2);  and  the  wake  with  initial  vertical 
momentum  Mq  is  assumed  to  be  a  dipole  of  constant  strength, 
initially  proportional  to  Mq,  with  a  moving  center.  The  use  of 
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such  point  sources  implies  that  hydrodynamic  disturbances  are 
not  expected  to  be  accurately  represented  at  short  distances 
from  the  source;  and,  inversely,  that  the  actual  source  size  can 
be  supposed  small  compared  to  a)  the  scale  over  which  the  density 
gradient  can  be  considered  constant,  and  b)  the  (vertical)  in¬ 
ternal  wave  length. 

As  mentioned  above,  the  ocean  environmental  conditions  are 
supposed  to  be  those  in  which  a  thin  thermocline  occurs.  This 
is  further  simplified,  in  the  model,  to  a  two  layer  fluid  with 
the  effects  of  the  actual  (assumed  weak)  density  gradient  below 

the  thermocline  is  assumed  to  be  incorporated  entirely  in  the 
source  strengths,  and,  for  moving  wakes,  in  geometrical  features 
of  the  trajectory.  Thus,  for  example,  the  strength  of  the 
quadrupole  representing  the  collapsing  wake  will  depend  on  the 
local  Vaisala  frequency  N,  and  on  features  of  its  geometry  and 
internal  turbulence,  but  the  disturbance  due  to  it  is  assumed 
to  be  equivalent  to  that  of  a  quadrupole  of  the  same  strength  and 
geometric  disposition  in  a  two  layer  fluid.  There  is  'supposed 
to  be  no  effect  on  source  strength  or  motion  of  the  deformation 
oi  the  density  interface,  which  will  generally  consist  of  a  local 
disturbance  and  interfacial  waves  radiated  away.  The  latter 
simulate  only  the  lowest,  "heaving"  mode  of  a  thin  thermocline; 
the  higher  internal  wave  modes  which  would  be  present  in  the 
actual  density  gradient  below  the  thermocline  are  simply  neglected. 

2 •  MATHEMATICAL  FRAMEWORK 

Across  the  density  interlace,  which  is  supposed  to  undergo 
small  displacements,  continuity  of  velocity  and  pressure  are 
assumed.  The  origin  of  coordinates  is  at  the  interface  with  the 
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axis  along  the  interface.  The  upper  layer,  of  density  px ,  has 
its  upper  boundary  at  y  =  h,  taken  to  be  a  rigid  wall;  the  bottom 
of  the  lower  layer,  of  density  p2,  is  at  z  =  -d.  The  source  of 
disturbance  is  at  a  distance,  ±f  from  the  interface.  The  medium 
is  supposed  to  be  at  rest  for  time  t  =  0.  Fourier  transform 
techniques  are  used  to  find  solutions  to  the  equations  of  motion 
and  boundary  conditions,  for  a  unit  "influence  function"  source 
at  x  =  §,  y  =  ±f  and  t  =  t.  Thus,  for  example,  the  velocity  at 
x,  h,  t  due  to  a  unit  source  at  (?,  -f,  t)  is 


u(x, t;  S,t)  = 


VST 


a 

I 


0s  sin  (s(x-?))ds 


CO 

/ 


+  I  sa9  sin  (a(t-T))  sin  (s(x-^))ds  [l] 

V  2tt 


where  5  and  H  are  the  delta  and  heaviside  unit  functions,  and 


2tt  cosh  ( 

Is |  (d-f)) 

! s !  sinh  ( 

fs(  (h+a)) 

g  j  s | e  sinh  ( 1 s 1 h )  sinh  (ls|d) 

sinh  ( | s I  (a+h) ) 


[2] 


For  a  source  function  having  prescribed  time  dependence  cp(T) 
the  velocity  U  at  the  wall  is 


U(  X,  h,  t )  =  f  u(  x,  t ;  5,t)  cp(  t  )  dT 


[3] 


HYDROWAUTICSj  Incorporated 


-5- 

The  results  for  (vertical)  dipoles  or  quadrupoles  are  obtained 
by  differentiating  with  respect  to  f.  Stationary  phase  approxi¬ 
mations  can  be  found  in  the  usual  way,  when  periodic  factors  can 
be  isolated  in  the  integrals  to  obtain  asymptotic  relations  for 
large  t.  For  other  cases  numerical  procedures  have  been  used. 

Thus  for  a  half-cycle  harmonic  source  at  (0,  -f)  whose 
strength  varies  with  time  >0  according  to  A  sin  yi }  and  is  0  for 
T  >  Tx  ,  where  y  =  2w/Tx  , 


Ti 


Ui  (x,  t) 


A  sin  y  t  u(x,t;  0,T)dT 


[4] 


o 


where,  for  t  >  Tx 


u(x,t;  0,t)  = 


Vf 


J 


a  9  sin  (a(t-T))  sin  sx  ds  [5] 


o 


The  asumptotic  formula  1  or  Uj  is  for  x  —  0,  x/t  —  fi  —  constant, 
(p 2  -  Pi}/ P3  =  e  «  1  and  large  t. 


Ui(jB,t) 


-2A 


2j  e 


-s  (f+h) 
ov  ' 


YOo 

— g - S'  sin 

Y -o 

o 


a  Tx  a  T 

-V  sin  (x„‘  +  -t +  ?) 


[6] 


xvhere  x(s)  -  /3s  -  cr(s),  and  SQ>  Oq  denote,  respectively,  the 
value  of  d3x/dss,  s  and  a  evaluated  at  sq,  the  root  of  da/ds  =  j8. 
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Similarly,  the  velocity  at  the  upper  boundary  due  to  a 
quadrupole  at  (0,  -f)  whose  strength  varies  with  time  as  Bte 
is 


Ua(x,t)  =  U2  H-  U2 


where 


-vt 

4Bte  x  [3(f+h)s  -  3x2 1 

f  ( f +h ) 2  +  x2]3 


[7] 


and 

OD 

u2  =  2B  f  e  s(f+h)  - - — ^  (cos  (sx-ct)  -  cos  ( sx  +  at)) 

oJ  ■  (Y2+a2)  2 

-  ya  [sin  (sx-at)  +  sin  (sx+at)]  ds  .  [81 

As  mentioned  above,  this  result  has  been  used  to  model  the 
early  stage  of  wake  collapse. 

For  the  final  stage  of  wake  collapse,  two  sources  and  two 
sinks  are  supposed  to  have  been  concentrated  at  (0,  -f)  at  t  =  0. 
As  time  increases,  the  two  sources  move  outwards,  in  opposite 
directions  with  constant  velocity  W,  while  their  strengths  vary 
as  M/ t  where  M  =  constant.  The  two  sinks  spread  out  to  form  a 
line  sink  of  uniform  strength  -  V/Wt2  joining  the  two  sources, 
(Reference^) .  For  this  case  the  velocity  induced  at  the  surface 
u3  =  U3  +  u3 
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where 


U-  = 


2M 


x-Wt 


+ 


x+Wt 


(f+h)2  +  (x+Wt)2 
(f+h)2  +  (x-Wt)2 


[9] 


is  a  secular  term,  and 


U,  = 


i  £ 

o  c 


s(f+h)  #  i  sin  sWt 

e  a  sm  s  x  sin  (o(t-T))  j  — — 


1 

COS  SWt  j  ds 

J 

[101 


3.  COMPARISON  WITH  EXPERIMENT 

Comparison  was  made  with  results  of  two-dimensional 
experiments  by  Wu  (Reference  2)  in  which  measurements  were  made 
of  interfacial  and  internal  waves  and  velocities  at  the  surface 
due  to  collapse  of  a  premixed  mass  of  fluid.  The  density  pro¬ 
file  in  these  experiments  involved  a  sharp ,  variable  density 
jump  above  a  constant  linear  density  gradient;  the  collapse 
occurred  at  a  fixed  location  relative  to  the  density  interface. 
Three  different  values  of  density  jump  were  used.  Wu's  observa¬ 
tions  indicated  that  for  the  largest  density  jump,  interfacial 
waves  clearly  predominated;  for  the  intermediate  value  of  density 
jump  both  interfacial  and  higher  mode  internal  waves  were  present 
and  for  the  smallest  jump,  the  higher  mode  internal  waves  pre¬ 
dominated.  It  was  felt  that  some  comparison  should  be  possible 
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with  his  results  for  the  case  of  the  largest  density  Jump.  For 
this  purpose  we  have  used  the  asymptotic  approximation  and 
assumed  a  quadrupole  source  to  represent  the  collapse  with 

■f* 

strength  varying  at  Bte“v  where  B  =  Nsa4,  a  =  radius  of  mixed 
region . 

From  the  slopes  of  the  rays  along  which  wave  energy  was 
concentrated,  some  estimates  can  be  made  of  the  frequency  dis¬ 
tribution  of  the  source  associated  with  the  wake  collapse.  For 
the  case  in  which  interfacial  waves  predominated,  the  value  of 
\j  was  estimated  in  this  way  to  be  v  ~  0.8N.  The  value  of  N  in 
the  experiments  was  about  0.7  sec-1.  We  then  find  (U2/  Y"geh) 

max 

0.03  to  compare  with  Wu's  observed  0.04  (corresponding  to 
Run  4  on  p .  33  of  Reference  2). 

The  agreement  is  probably  somewhat  fortuitous  owing  to  the 
fact  that  the  quadrupole  strength  is  only  approximate.  However, 
the  maximum  U2/ ^fgeh  does  not  seem  sensitive  to  exact  values  of 
the  parameters  involved. 

If  the  same  procedure  is  followed  for  Wu's  case  2,  with  an 
intermediate  density  jump,  the  agreement  is  poorer.  Wu  found  a 
maximum  U2/ygeh  about  four  times  larger,  but  above  expression 
gives  a  smaller  value  than  for  the  previous  case.  This  dis¬ 
crepancy  can  be  ascribed  to  the  presence  of  the  higher  mode  in¬ 
ternal  waves  observed  by  Wu,  and  which  are  not  treated  by  this 
theory. 

3. a  Comparison  with  Miles'  Theory  for  Sharp  Thermocline  Case 

Comparison  can  be  made  with  the  theoretical  results  of  the 
three-dimensional  calculations  by  J  Miles  (Reference  4)  for  the 
case-  of  a  thin  thermocline.  In  a  plane  perpendicular  to  the 
direction  of  motion,  the  two-dimensional  model  of  the  hull  dis- 
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turbance  corresponds  to  flow  induced  by  a  time-dependent  monopole 

of  strength:  A  sin  yt  where  A  ~  a  2  V/L  in  which  a  =  radius 

o  o 

and  L  the  length  of  the  (assumed  axisymmetric )  hull,  V  the 
(constant)  forward  speed,  and  y  =  2tt/T,  with  T  =  L/V  . 

For  the  lowest  internal  wave  mode  of  a  thin  thermociine, 

the  asymptotic  two-dimensional  result  for  the  surface  strain  Ix 

due  to  the  hull  at  depth  f  +  h  is  given  by  Equation  [6]  above. 

Omitting  the  periodic  dependence  on  t,  and  following  the  same 

approximation  procedure  as  Miles',  the  maximum  Ix  is  found  to 

Q  1 

occur,  approximately,  at  sq  ~  J  for  which  (replacing  t  by 

I/V 

_  i +  -9/4 

|lilmax  "  °-4  a0’L  {-  V  (eO  (f+h)  I’ll] 


Except  for  notation.  Equation  [ll]  is  the  same  as  Miles'  result 
for  the  hull  disturbance  (Reference  4). 

The  asymptotic  expression  for  the  surface  strain  induced 
by  a  collapsing  wake,  modeled  by  a  two-dimensional  quadrupole 
in  a  transverse  plane,  with  the  time  dependence  Bte~vt  is 


2B 


Yx- 


-so(f+h) 


O 


lx"'2  (yz+o  2)2 
o  o' 


y2-o  2 

2~~  COS  (xot  +  ?) 


+  ya  sin  (xt  +  J) 


o 


o 


[12] 


11 


o  4  f+h  * 


With 


the  maximum  I2  is  found  to  occur  at  about  s  = 

o 

B  ~  W2a4 .  where  N2  is  the  Vaisala  frequency  at  depth  f+h  and 
a  is  the  effective  mixed  wake  radius  at  collapse: 
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/V  \s  -11/4  -3/4 

!  1=  !max  '  Y  (f+h)  (g€)  r  13] 

Except  for  notation  and  the  numerical  factor  in  front 
being  an  order  of  magnitude  larger,  Equation  j[l3j  is  the  same 
as  Miles’  result  for  the  collapsing  wake.  Despite  the  uncer¬ 
tainty  in  source  functions,  this  discrepancy  may  reflect  the 
fact  that  the  two-dimensional  approximation  tends  to  overestimate 
and  on  the  other.  Miles'  use  of  a  small-body  type  approximation 
tends  to  under-estimate  the  surface  disturbances  due  to  wake 
collapse . 

4.  APPLICATIONS  OF  MODEL 

This  model  has  been  employed  to  investigate  effects  due  to 
a  submarine  in  the  upper  mixed  layer,  above  a  thin  thermocline, 
but  running  out  of  trim,  and  "heavy".  Its  wake  will  then  have 
a  negative  momentum  in  the  vertical  direction,  moving  down  to¬ 
ward  the  thermocline.  The  wake  will  entrain  water  on  its  way 
down,  and  slow  down;  however,  we  have  assumed,  for  simplicity, 
a  constant  vertical  dipole  moment  and  constant  velocity  for  the 
wake  until  it  reaches  the  thermocline.  Actually,  the  velocity 
is  inversely  proportional  to  the  distance  traveled  by  the  wake, 
divided  by  the  initial  effective  wake  radius  (Reference  5);  under 
the  geometrical  circumstances  assumed  the  velocity  probably 
would  not  vary  by  as  much  as  50  percent.  When  the  wake  encounters 
the  thermocline,  it  is  assumed  the  downward  motion  is  arrested, 
simulated  by  the  addition  of  a  negative  dipole  moment,  after 
which  the  dipole  moment  is  assumed  to  undergo  a  cycle  of  damped 
oscillation  to  zero,  similar  to  phenomena  noted  in  Reference  5. 
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Before  all  this,  one  cycle  of  a  two-dimensional  monopole  is 
assumed  to  occur;  representing  the  effect  of  the  hull.  The 
geometrical  arrangement  and  time-history  of  these  sources  is 
depicted  in  Figure  1.  ,  The  hydrodynamic  source  parameters  were: 
representing  the  hull,  a  monopole  source  =  A  sin  yt  with 

A  ~  a02VL  and  Y  ~  V0/L’  with  Vo  ~  6  knots;  radius  aQ  =  103  cm, 

somewhat  larger  than  usual,  and  length  to  radius  ratio  L/a  ~  IB 
*  o 

The  depth  is  20  meters  above  the  thermocline,  which  is  at  60 

meters.  Across  the  thermocline  the  relative  change  in  density 
&p/p  =  3.10-4.  The  wake  dipole  moment  has  the  initial  value 
'i  da”o  3  corresPonding  to  a  few  tons  "heavy".  is  assumed 

to  remain  constant,  with  its  center  moving  downward  at  a  constant 
speed  5  cm/sec,  until  within  aQ  cm  of  the  interface,  when  Dx  is 
switched  off  by  adding  D2  =  -D1 .  D3  is  then  switched  on,  repre¬ 
senting  the  overshoot  and  fail-back  of  the  wake  by  the  damped 
oscillatory  function  Dx  (cos  at)e~at  with  a  ~  3/T,  T  =  (ge/h)"^  ~ 
140  sec  is  the  daisaia11  period  of  the  interface  itself. 

Figure  2  shows  the  behavior  of  the  velocity  at  the  wall 
for  each  of  the  sources  of  Figure  1^  separately ^  and  for  the 
superposition  of  all  sources  together.,  for  two  multiples  of 
interface  Vaisala  periods  2.8T  and  5*6t^  corresponding  roughly 
to  lOO  and  800  sec.  Evidently  the  sum  of  all  the  contributions 
to  the  velocity  at  the  wall  exhibits  effects  of  some  destructive 
interference  being  smaller  than  the  contributions  for  most  of 
the  individual  sources.  None  of  the  individual  sources  is 
negligible  compared  to  the  rest.  Undoubtedly  the  destructive 
interference  reflects  the  particular  assumptions  involved  in 
this  calculation.  The  time  phasing  assumed  depends  on  the 
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assumed  trajectory  of  the  downward  moving  dipole.  Actually,  It 
would  be  expected  that  the  wake  trajectory  will  differ  from  that 
assumed,  due  to  the  entrainment,  and  would  slow  down  as  it  moved 
toward  the  thermocline  (Reference  5)-  Also  the  effective  radius 
of  the  wake  will  grow,  so  that  it  will  reach  the  thermocline 
somewhat  sooner  than  assumed;  the  dipole  effects,  therefore,  may 
be  somewhat  exaggerated.  Calculations  taking  into  account  these 
adjustments  have  not  been  carried  out. 

Corresponding  to  the  situation  in  which  a  submarine  is  run¬ 
ning  slightly  "light"  below  a  thin  thermocline.  Figure  3  shows 
the  time  sequence  assumed  for  wake  with  upward  momentum  and  wake 
collapse.  The  wake  is  assumed  to  have  a  damped  upwards  trajec¬ 
tory  oscillating  about  the  same  level  at  which  it  collapses 
''Reference  6),  ten  meters  below  the  density  interface,  which  is 
at  a  depth  of  50  meters.  Shallow  conditions  are  assumed,  with 
the  bottom  at  150  meters.  The  source  parameters  correspond  to 
a  wake  radius  a  =  103  cm,  and  a  maximum  upwards  velocity  of 
3  cm/sec,  corresponding  to  a  dipole  moment  D  =  3.10s  (jos  at)eat, 
with  a  =  1/4T,  where  T  =  the  interface  Vaisala  period,  again 
about  140  sec.  The  quadrupole  strength  is  taken  to  be  Q  =  Bte_v 
with  B  -g-I^a^4,  where  now  N  is  the  Vaisala  frequency  of  the 
actual  density  gradient  below  the  thermocline  and  27t/N  ~  4t, 
y  =  1/4T.  The  velocities  at  the  surface  due  to  the  dipole  and 
quadrupole  separately  and  superposed  are  shown  in  Figure  4  as  a 
function  of  distance  at  t  =  2ttT  ~  8l6  sec.  In  this  case  the 
quadrupole  appears  to  have  a  larger  effect,  by  about  a  factor  2, 
and  the  individual  contributions  reinforce  each  other. 


HYDRONAUTICS,  Incorporated 


-13- 

The  relative  separation  in  time  of  these  sources  shown  in 
Figure  3  may  be  larger  than  actually  is  the  case.  Tulin  has  ob¬ 
served  that  the  collapse  often  begins  at  the  same  instant  as  the 
wake  begins  to  fall  (Reference  6).  If  this  would  be  the  case, 
the  quadrupole  effects  would  be  larger.  Laboratory  experiments 
reveal  a  complex  behavior  of  the  terminal  phase  of  trajectories 
of  turbulent  vortex  pairs,  simulating  wakes  with  vertical  mo¬ 
mentum  in  a  density  gradient;  and  the  final  wake  collapse ,  in 
such  cases,  appears  to  be  a  more  complicated  function  than  shown 
in  Figure  3  (Reference  5)*  Calculations  with  different  values 
of  relative  source  phasing  have  not  been  carried  out. 

It  may  be  worth  noting  also  that  the  results  shown  in 
Figures  2  and  4  are  derived  from  numerical  evaluation  of  integral 
expressions  for  the  full  disturbance,  and  not  only  the  asymptotic 
approximations . 

The  model  has  been  employed  also  to  explore  the  effect  of 
different  ratios  of  bottom  depths  to  thermocline  depths.  The 
bottom  effects  on  disturbances  due  to  the  hull  are  generally 
small  except  for  quite  shallow  depths.  Such  effects  for  the  quad¬ 
rupole,  for  the  same  assumed  conditions,  are  generally  smaller. 

Another  problem  to  which  this  model  has  been  applied  is  that 
of  the  "efficiency"  of  internal  wave  production  by  a  collapsing 
wake.  Efficiency  is  here  defined  as  the  ratio  of  the  energy  of 
the  lowest  mode  internal  wave,  long  after  collapse,  to  the 
initial  potential  energy  of  the  collapsing  mass  of  fluid,  assumed 
fully  mixed  and  initially  quiescent.  For  such  a  fixed  time,  the 
energy  of  the  lowest  mode  internal  wave  was  estimated  as 
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oo 

2  j  Ap  y2gdx;  where  Ap  is  the  density  jump  across  the  thin 
J 
o 

thermocline  and  y  is  the  thermocline  displacement.  The  time 
dependence  of  the  collapsing  wake  was  again  taken  ~Bte_V^:  where 
B  was  taken  to  be  N2a4  where  N  is  local  Vaisala  frequency  and  a 
the  effective  wake  radius.  Total  potential  energy  of  the  wake 
was  taken  to  be  ir/4  N2a4  .  The  result  for  Wu's  experiments  was 
about  6  percent,  which  seemed  surprisingly  small.  A  check  on 
this  was  made  using  Wu's  experimental  values  for  interface  dis¬ 
placement  to  evaluate  the  above  integral,  and  gave  nearly  the 
same  result.  Lighthill  had  remarked  that  the  efficiency  of  low¬ 
est  mode  generation  by  a  collapsing  wake  would  be  small  (Refer¬ 
ence  7 >  see  also  Reference  9) • 

Some  exploration  was  also  conducted  of  effects  of  different 
time-dependencies  of  the  quadrupole  source.  First  by  using  the 
same  functional  form,  i.e.,  strength  proportional  to  t  e_vt,  and 
calculating  surface  strains  and  currents  for  fixed  source  geome¬ 
try  and  values  of  density  jump  across  the  thermocline,  for  a 
range  of  values  of  v.  Some  appreciable  variation  of  the 
was  found.  Next  calculations  were  made  for  a  more  general 
functional  form  tne"vt,  where  n  was  allowed  to  be  <T  and  >1.  This 
investigation  could  not  be  completed. 

It  may  be  noted  that  the  secular,  local  disturbance  velocity 
induced  by  a  quadrupole  of  strength  Bte"Mt  located  at  depth  f+h 
below  a  thermocline  if  depth  h  is  given  by  U2  from  Equation  [j] 
above.  Near  x  =  0,  the  surface  strain  due  to  wake  collapse, 
modeled  by  this  source,  approaches,  for  long  times 
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I2-  -  12B/vs(f+h)4 

or,  with  B  ~  N2a4,  I2-*  12N2a4/v2  (f+h)4  ;  with  v  ~  N,  if  the 
maximum  effective  wake  radius  a  is  of  the  order  of  4  times  the 
original  wake  radius  a  ,  Reference  (8),  quite  appreciable  values 
of  I2  can  result.  A  more  accurate  model  of  final  wake  collapse 
at  large  times  would  involve  the  source-sink  combination  dis¬ 
cussed  above.  Equation  [9]j  this  has  not  been  carried  out  but  it 
is  believed  that  doing  so  would  increase  the  value  of  I2  above. 

5 .  CONCLUDING  REMARKS 

It  should  perhaps  be  emphasized  again  that  the  results 
given  above  involve  only  the  lowest  mode  internal  wave  of  a  thin 
thermocline.  The  fact  that  it  appears  to  be  the  most  persistent 
in  usual  ocean  conditions  lends  importance  to  this  mode.  However, 
for  thermoclines  of  thickness  b  >  a  ,  and  submarine  depths  in  the 
thermocline,  the  second  and  higher  modes  are  likely  to  be  more 
important.  It  would  be  desirable  to  extend  this  type  of  calcu¬ 
lation  to  include  these  higher  modes.  Because  of  the  differences 
in  speeds  of  internal  wave  modes  and  of  their  dispersive  char¬ 
acteristics,  the  results  can  be  expected  to  be  correspondingly  more 
complicated. 

Surface  effects  associated  with  the  lowest  mode  of  a  thin 
thermocline  have  a  dependence  on  depth  dominated  by  the  ex¬ 
ponential  in  Equation  [61  which  can  be  regarded  as  the  product 
of  two  factors:  one  a  function  of  depth  below  the  thermocline 

exp(-s  f)  where  s  is  the  wave  number  of  the  lowest  mode,  and 
* K  o-l  ’  ox  3 
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the  other,  exp(-SQ^h),  of  depth  above  the  thermocline.  The  depth 
dependence  of  higher  modes  for  a  thick  thermocline  with  a  more 
or  less  well-defined  upper  edge  at  depth  h  will  characteristically 
have  an  exponential  dependence  on  h,  exp(-sQ  h)  where  s  is  trie 
(larger)  wave  length  of  the  nth  mode,  but  a  more  or  less^periodic 
dependence  on  depth  below  the  thermocline.  Often  the  upper  edge 
of  a  thick  thermocline  is  not  sharply  defined,  however,  and  h  is 
then  somewhat  different  for  the  various  modes. 

Speed  dependence  for  the  hull  effects  may  be  expected  to 
increase  at  low  speeds  because  of  the  'resonance"  type  denominator 
in  Equation  [6].  With  a  similar  source,  the  same  type  of  denom¬ 
inator  may  be  expected  to  appear  for  every  mode,  with  appropriately 
changed  values  of  .  At  higher  speeds,  the  hull  contribution 
decreases  as  in  Equation  [ll].  The  dipole  strength  will  depend 
on  speed  through  the  vertical  velocity ^  which  is  proportional  to 
the  forward  speedy  for  the  same  angle  of  attack  of  the  sailplanes 
or  hull  (Reference  1) . 

The  quad rup ole  contribution  increases  with  speed  through  the 
dependence  on  V  explicit  in  Equation  [13]  and  implicitly  through 
dependence  of  effective  wake  size  on  speed.  In  a  linear  density 
gradient^,  rough  upper  bounds  to  this  wake  size  can  be  obtained  by 
assuming  all  the  work  of  turbulent  expansion  of  the  wake  from  a 

o 

to  a  goes  into  increasing  potential  energy.  This  gives  (a/a  )4  =* 

V  2  v  o' 

-j  1  o 

1  +  TjT  2  which  seems  to  agree  with  most  of  the  available 
o 

laboratory  data.  However,  for  our  model  this  can  only  apply  for 
f  ~  a.  For  smaller  f,  the  wake  encounters  the  thermocline  during 
its  growth,  and  the  expansion  beyond  the  thermocline  must  take 
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into  account  the  nonlinearity  of  the  overall  density  gradient. 

For  higher  modes  in  a  thick  thermocline,  this  type  of  speed  de¬ 
pendence  will  hold  for  the  modes  having  vertical  wavelengths 
much  larger  than  the  size  of  the  wake;  this  depends  on  the  thick¬ 
ness  of  the'  thermocline.  The  modes  with  shorter  wavelengths 
should  not  have  such  strong  speed  dependence  (Reference  9) • 

Because  of  the  propensity  for  internal  waves  to  have  non¬ 
linear  interactions ,  superposition  of  the  internal  wave  distur¬ 
bances,  as  assumed  above,  may  exaggerate  the  surface  effects  so 
estimated . 

The  hydrodynamic  source  functions  have  been  assumed  to  be 
"decoupled"  from  the  internal  wave  field.  Some  basis  for  this 
assumption  can  be  found  in  the  fact  that  the  trajectory  of  tur¬ 
bulent  vortex  pairs,  and  at  least  some  aspects  of  wake  collapse 
phenomena  can  be  estimated  without  having  to  deal  with  the 
coupled  phenomena  (References  2  and  5)-  An  attempt  was  begun  to 
obtain  effective  (uncoupled)  source  strengths  by  systematic 
comparison  of  experimental  data  on  interface  deformations  (Ref¬ 
erences  2  and  5)  with  calculations  involving  the  functional 
forms  of  sources  assumed  above,  in  order  to  obtain  "best  fit" 
parameters  in  these  functions.  Unfortunately,  this  work  could 
not  be  completed. 

It  is  felt  nonetheless  that  the  results  indicate  (1)  that 
the  effects  of  wake  vertical  momentum  should  be  taken  into 
consideration  in  the  estimation  of  surface  effects  due  to  passage 
of  a  submerged  submarine  in  a  stratified  ocean;  and  (2)  that  the 
superposition  of  the  complex  pattern  of  disturbances  associated 
with  actual  submarine  passage  can  give  rise  to  significant 
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variation  in  the  magnitude  of  surface  effects.  It  is  recognized 
that  these  results  are  rather  fragmentary;  it  is  hoped,  however, 
that  they  are  sufficiently  suggestive  to  stimulate  further  work 
with  more  accurate  models,  over  wider  ranges  of  operational  and 
environmental  parameters,  and  taking  into  account  to  the  extent 
possible  the  coupling  of  the  disturbance  field  to  the  dynamics 
of  the  wake . 


HYDRONAUTICS,  Incorporated 


-19- 


REFERENCES 

1.  Tulin,  M.  P.,  in  Proceedings  of  the  Surface  Effects 
Symposium,  1964.  (SECRET) 

2.  Wu,  J.,  "Flow  Phenomena  Caused  by  the  Collapse  of  a  Mixed 
Region  in  a  Density-Stratified  Medium, "  HYDRONAUTICS, 
Incorporated  Technical  Report  231-11,  1966 .  (Unclassified) 

3.  Tulin,  M.  P.,  "Hydrodynamic  Effects  and  Detection  (U)," 
HYDRONAUTICS,  Incorporated  Technical  Report  231-2,  1963* 
(SECRET) 

4.  Miles,  J.,  Jason  Report  S334,  1968.  (SECRET) 

5.  Tulin,  M.  P.,  and  Schwarz,  J.,  "The  Motion  of  Turbulent 
Vortex  Pairs  in  Homogeneous  and  Density-Stratified  Medium, " 
HYDRONAUTICS,  Incorporated  Technical  Report  231-15 *  1966 . 
(Unclassified) 

6.  Tulin,  M.  p..  Proceedings  of  the  Non-Acoustic  Surface 
Effects  Symposium,  1968.  (SECRET) 

7.  high thill,  M.  J.,  SPOSS  Report,  RRE  Malvern,  1964.  (SECRET) 

8.  Van  de  Watering,  W.P.M.,  Tulin,  M.  P.,  and  Wu,  J.,  "Experi¬ 
ments  on  Turbulent  Wakes  with  Stable  Density  Gradient, " 
HYDRONAUTICS,  Incorporated  Technical  Report  231-24,  1971* 

( Unclassified) 

9.  Schooley,  A.  H.,  and  Hughes,  B.,  "An  Experimental  and 
Theoretical  Study  of  Internal  Waves,  Generated  by  the 
Collapse  of  a  Two-Dimensional  Mixed  Region  in  a  Density 
Gradient,"  Journal  of  Fluid  Mechanics,  Vol.  51*  1972. 
(Unclassified) 


HYDRONAUTICS,  INCORPORATED 


TIME ,  seconds 

I  1 

0  100  200  400  600  800  1000  1200  1400 

I _ l _ i _ i _ i _ i_ _ l _ i _ i _ | 


+ 


TIME,  T  (  Interface  Vaisala  Periods  ) 


9 


10 


FIGURE  1  -  TIME  HISTORY  OF  HYDRODYNAMIC  SOURCE  CONTRIBUTIONS  SIMULATING 
SUBMARINE  RUNNING  HEAVY  ABOVE  THERMOCLINE. 
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SURFACE  DUE  TO  HYDRODYNAMIC  SOURCES  SHOWN  IN 
FIGURE  1  AT  2.8  AND  5.6  INTERFACE  VAISALA  PERIODS  T. 
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FIGURE  3  -  TIME  HISTORY  OF  HYDRODYNAMIC  SOURCES  SIMULATING  WAKE  OF 
SUBMARINE  RUNNING  LIGHT  JUST  BELOW  THERMOCLINE. 
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FIGURE  4  -  VELOCITIES  AT  SURFACE  DUE  TO  SOURCES  IN  FIGURE  3 
SIMULATING  WAKE  OF  SUBMARINE  SLIGHTLY  LIGHT 
BELOW  THERMOCLINE 
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